Lagrangian tracers on a surface flow: the role of time correlations 
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Finite time correlations of the velocity in a surface flow are found to be important for the formation 
of clusters of Lagrangian tracers. The degree of clustering characterized by the Lyapunov spectrum 
of the flow is numerically shown to be in qualitative agreement with the predictions for the white- 
in-time compressible Kraichnan flow, but to deviate quantitatively. For intermediate values of 
compressibility the clustering is surprisingly weakened by time correlations. 
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Inhomogeneous distribution of particles advected in a 
turbulent flow is a generic consequence of compressibil- 
ity. This can be obtained in two situations. The first 
possibihty is that the advecting flow is compressible it- 
self and that the particles follow the streamlines [1]. The 
other possibility is that the particles do not follow the 
streamlines because of inertia [2, 3, 4, 5, 6, 7, 8] or lift 
[9], and that the effective velocity field is compressible. 
Such situations are relevant for the formation of clouds 
[10] or for the advection of bubbles in turbulent flows, 
e.g. for breaking waves on the ocean surface [11, 12]. We 
focus here on the first possibility. 

While the dominant tendency of incompressible flows 
is to separate particle trajectories, a compressible com- 
ponent is responsible for particle trapping in contract- 
ing regions for long times. The Eulerian compressibil- 
ity of a flow is measured by the dimensionless ratio 
C = {{diUif') / {{diUj)^) . It takes values between (in- 
compressible fiow) and 1 (potential fiow). While there 
can be no clustering without compressible effects, the 
compressibility ratio C is insufficient to determine the fi- 
nal distribution completely [13]: the behavior depends 
also on the spatial roughness, the dimensionality, and, as 
we will demonstrate here, on the time correlations in the 
fiow. 

A convenient characterization of the final distribution 
uses the Lagrangian Lyapunov spectrum. Dynamical 
systems theory shows that the asymptotic clusters are 
smooth along the unstable directions of positive Lya- 
punov exponents and fractal along the stable directions. 
In d dimensions, the sum of Lyapunov exponents 
vanishes in the incompressible case and is negative for a 
compressible fiow. A measure for the distribution of the 
final clusters is given by the Lyapunov dimension [14], 



Dl = K 
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where K is the maximal integer such that X^i^i — ^■ 
Obviously, for an incompressible fiow K — d and also 
Dl = d: the particle distribution fills the entire volume. 
As compressibility increases, the sum of Lyapunov ex- 
ponents will become negative and K will drop below d. 



If the largest Lyapunov exponent Ai becomes negative, 
the final cluster will not have any smooth directions any- 
more and the particles will cluster in a point-like fractal. 
One hence distinguishes a regime of strong compressibil- 
ity Ai < with Dl = 0, and one of weak compressibility 
Ai > 0. 

There are no general results on the spectrum of La- 
grangian Lyapunov exponents in turbulent flows. For 
the case of incompressible, isotropic, three dimensional 
turbulence the numerical observation is that A2 ~ 
l/4Ai [15]. In the limit of a compressible Kraich- 
nan fiow, which is a synthetic, white-in-time, and 
Gaussian distributed, the spectrum is given by Aj = 
Ci [d{d ~ 2j + 1)- 2C{d +{d- 2)j)] where j = 1, d, 
Ci is an inverse time proportional to the Lagrangian 
strain and the resulting Lyapunov dimension is a decreas- 
ing function of C [13]. Moreover, when C > Cg = d/A all 
the Lyapunov exponents become negative and one has 
the strong compressible regime. 

The trapping effects are believed to be enhanced in 
spatially rough Kraichnan flows [13]. Although a recent 
attempt has been devoted to re-introduce the flnite time 
correlation for the synthetic rough Kraichnan flows [16] a 
general theoretical framework studying the effects of the 
time correlation is still lacking. 

Here, we demonstrate numerically the crucial role of 
finite time correlations in the Lagrangian statistics of 
particles transported by the Navier-Stokes fiow which is 
established in a free-slip surface. As a consequence of 
boundary conditions the two-dimensional fiow displays 
an effective compressibility. By decomposing the surface 
velocity field in its irrotational and potential components, 
we change the effective value of compressibility of the sur- 
face fiow. This enables us to study the interplay of the 
two components, time correlations and compressibility, 
within one system. 

We want to remark that our investigation does not en- 
counter the effects of fiuid density variations which ap- 
pear in a compressible turbulent fiow. 

The tracers are advected by the free surface flow on top 
of three-dimensional (3c?) incompressible Navier-Stokes 
turbulence at Re\ ~ 145. The volume with an aspect 
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FIG. 1: Compressibility C as a function of depth zjLz rang- 
ing from the free surface {z/ = 0) to the middle of the tank 
{z/Lz = 0.5). The figure illustrates also that for z/L > 0.1 
the volume turbulence gets close to the homogeneous isotropic 
case which is indicated by the vertical dotted line. The hor- 
izontal dotted line indicates the C value of 1/6 which would 
follow for perfect isotropy in the bulk. 



ratio of 27r : 27r : 1 is resolved by 256 x 256 x 65 grid 
points. The full Sd equations are integrated by a stan- 
dard pseudo-spectral method [17]. The lateral bound- 
ary conditions for x and y are periodic and in the ver- 
tical direction we apply free-slip boundary conditions, 
Mz = and dzUx — d^Uy ~ 0. The Lagrangian trac- 
ers are advected by the surface flow spanned by the two 
components Ux{x, y, z = 0, t) and Uy{x, y,z = 0, t) which 
is compressible since dxU^ -\- dyUy = —dzUz ^ [17]. 
Tracer velocity between the grid mesh is calculated by 
bi-linear interpolation. We store 18,000 configurations 
of the 2d surface velocity field at the same resolution in 
time that is used for the time advancement of the full 3d 
Navier-Stokes equations. On the basis of these configu- 
rations, Lagrangian trajectories and the Lyapunov spec- 
trum (Ai, A2) are computed. 

Isotropy can be used to show that the compressibility 
ratio for a two-dimensional flow {ux{x,y),Uy(x,y)) can 
be rewritten as 



C 



{{dxUx + dyUy)' 

2{{dxUx - dyUy) 



(2) 



Numerical and experimental investigations of the surface 
flow at Taylor Reynolds numbers Re\ ^ 10^ have shown 
that C ~ 0.5 [18, 19]. By (2), this is equivalent to a 
vanishing value of the cross correlation (dxUxdyUy), as 
confirmed in numerical studies [20]. Figure 1 shows the 
profile of C for all planes of the simulation box from z = 
down to z = L/2. It is interesting to note that the value 
of compressibility measured in the bulk is very close to 
the value 1/6 which is expected for a two-dimensional 
cut of a three-dimensional homogeneous, isotropic flow. 
The relatively large value of the compressibility ratio on 
the free surface is a consequence of the free-slip boundary 
conditions. 



The evolution of an initially uniform particle distribu- 
tion in the time-correlated surface flow with C = 0.45 is 
shown at two instants in Figs. 2a and 2c, respectively. 
The pronounced formation of particle clusters was ob- 
served recently in experiments [18] and in simulations 
[19]. 
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FIG. 2: Visualization of the particle clustering for a uniform 
initial distribution of 36,000 tracers on the free surface z = 
with C = 0.45. The panels (a) and (c) show the tracer 
distributions in time correlated surface flow at 0.2r and 3. IT. 
The panels (&) and (d) correspond to the tracer distributions 
in the time de-correlated fiow at 0.2T and 5. IT. The time 
scale T = L/5u{L) is the large scale eddy turnover time. 

Although having a compressibility close to the critical 
value of C2 = 1/2 for the two-dimensional Kraichnan flow 
we do not observe point-like structures; particles accu- 
mulate on a network of narrow ridges with bigger empty 
voids in between. Figure 3 reports the results of the nu- 
merical computation of the Lagrangian Lyapunov spec- 
trum. The very existence of the ridge structures pertains 
to the existence of a finite positive Lyapunov exponent 
Ai ^ 0.3. This is in line with the computed Lyapunov 
dimension, Dl — 1.15 which is plotted in the inset of 
Fig. 3. 

We now compare the Lyapunov dimensions of the time- 
correlated flow with the the corresponding de-correlated 
flow. That would quantify the effects of time correla- 
tion effects on clustering properties. In order to do so, 
we generate the decorrelated velocity field from the flow 
at hand by re-shuffling the temporal sequence of surface 
flow snapshots. Thanks to the periodic boundary con- 
ditions, each snapshot is further de-correlated by phase 
randomization in space. 

Figures 2b and 2d show the evolution of the density of 
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FIG. 3: Convergence of Lyapunov exponents (Ai > and 
A2 < 0) for the Lagrangian trajectories on the free surface 
computed for about three eddy turnover times T. In the 
inset, the Lyapunov dimension Dl = 1 + A1/IA2I is shown. 



tracers in the re-shuffled time de-correlated flow. We find 
the clusters concentrated on more point-like sets and the 
ridges are depleted evidently. Therefore one expects that 
the Lyapunov dimension is decreased when the temporal 
correlation of the velocity is diminished. Lagrangian tra- 
jectories in the randomized flow are more erratic: in order 
to achieve statistical convergence, trajectories are now in- 
tegrated for about 5 large scale eddy turnover times. The 
computation of Lyapunov spectrum and dimension for 
the randomized flow is shown in Figure 4. The fact that 
both of the Lyapunov exponents are changed with respect 
to Fig. 3 is not relevant, as the characteristic Lagrangian 
time is changed by time decorrelation procedure. The 
important point is that the Lyapunov dimension is now 
Dl — 1.05 which is less than the Lyapunov dimension in 
the time correlated flow. This is quite surprising as one 
intuitively expects that the particle trapping is enhanced 
when moving in a time-correlated flow because the local 
flow will persist for some time. 

In the case of two-dimensional compressible Kraichnan 
ensemble the pairing symmetry in the Lyapunov spec- 
trum (A, —A) is broken and the Lyapunov dimension is 
then simply given by 



1 + 2C 



(3) 



According to (3), the transition to the strong compress- 
ibility regime takes place at Cg = 1/2, very close to 
the numerically observed compressibility at the surface 
(see Fig. 1). The prediction of (3) for C ~ 0.45 renders 
Dl ~ 1.06 which is in very good agreement with the 
result of randomized surface flow. The compressibility, 
seen along Lagrangian trajectories [4], is different from 
C as given by (2). However explicit theoretical results 
inferring to the effects of velocity time correlations on 




FIG. 4: Convergence of Lyapunov exponents for the La- 
grangian trajectories advected by the random-time, random- 
phase scrambled velocity field. The computation is done for 
5T. In the inset, wc show the Lyapunov dimension Dl to- 
gether with the theoretical prediction Dl = 1.05 (dashed 
line). 



effective Lagrangian compressibility or Lyapunov spec- 
trum are still lacking although first attempts have been 
made [16]. Moreover, it is not clear that the theoreti- 
cal description of the time correlated flows will enjoy the 
same level of generality. 

In order to gain more insight on the role of the time 
correlations we repeated the Lagrangian analysis for dif- 
ferent degrees of compressibility. To this end, wc have 
decomposed the 2d velocity snapshots into the stream 
function il){x,y) and the potential function (f){x,y). A 
new velocity field is reconstructed then to 

u{x,y) = V2 (y-^ipix, y) cos 7 -|- V(j){x, y) sin 7) (4) 

where V"^ = {dy,—dx)- The compressibility ratio C 
(cf. (2)) for the decomposed flows is a monotonic function 
of the parameter 7 G [G,7r/2] and for 7 = 7r/4 the origi- 
nal surface flow is reproduced. 7 = corresponds to the 
incompressible case and 7 = 7r/2 gives a purely potential 
flow; so we are able to cover the whole range of compress- 
ibility degrees. Of course, dimensional quantities such as 
the characteristic Lagrangian time of the velocity field 
may depend on the value of 7. On the other hand, di- 
mensionless quantities, such as Dl are independent on a 
global rescaling of the velocity and thus can be expected 
to depend mainly on the degree of compressibility. 

We have performed two sets of numerical simulations of 
Lagrangian trajectories for different values of 7. For each 
value, we have computed the mean compressibility and 
the Lyapunov dimension, both for the time-correlated 
and randomized flows. The result is summarized in Fig. 5 
where for comparison we plot the relation (3) for the S- 
correlated Kraichnan case. 
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the tracer distribution is smaller than the dimension in 
the time de-correlated flows. Yet as the compressibility 
gets higher this regime crosses over to a more interesting 
one where the time correlations increase the Lyapunov 
dimension. The point of phase transition to strong com- 
pressible regime is at a higher value, C ~ 0.75. The 
present paper is a first attempt to relate aspects of the 
solvable compressible Kraichnan model to the turbulent 
Navier-Stokes compressible flows. This suggests further 
investigations on the interrelated roles of the time corre- 
lations and the compressibility on the dispersion proper- 
ties. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. The computations were carried out at 
the John von Neumann Institute for Computing at the 
Forschungszentrum Jiilich and we are grateful for sup- 
port. G. B. acknowledges hospitality of the Philipps 
Universitat Marburg. 



FIG. 5: Lyapunov dimension Dl computed for different val- 
ues of the compressibility C in accordance with (4). Data 
points for time-correlated flow are indicated by squares. The 
open circles are for the de-correlated flow. The line indicates 
the behavior for the compressible Kraichnan flow (cf. (3)). 
The gray shading in the background illustrates the border 
between weak and strong compressible regime for the time- 
correlated flow. We have also indicated the data points in the 
plane that correspond with the tracer distributions as given 
in Fig. 2. 



It is remarkable that the qualitative effect of time cor- 
relations on tracers distribution depends on the level of 
compressibility. For small values C < 0.3, time corre- 
lations enhance the effects of compressibility, i.e. Dl 
smaller than the predicted value by (3). The situation 
for larger values of C is more surprising as the Lyapunov 
dimension in correlated flow is larger than the value given 
by (3), i.e. the effects of compressibility are depleted 
with respect to the i5-correlated case. The later situation 
is contrary to intuition according to which finite time 
correlations amplify the effects of compressibility. Con- 
sequently the transition point to the strong compressible 
regime, where Ai < 0, moves up to the C ~ 0.75 larger 
than the prediction of (3) (cf. gray region in Fig. 5). 
This is again the point where the ratio gets zero for 
the correlated flow. 

In summary, we have studied the effects of time corre- 
lations in the clustering of Lagrangian tracers advected 
on a free surface flow by direct numerical simulation 
of NS equations. Appealing to a numerical decompo- 
sition method we generate 2d flows with variable degrees 
of compressibilities. The Lyapunov dimension decreases 
with the compressibility degree. Our simulations indicate 
that the effect of time correlation can go both ways. At 
low values of compressibility the Lyapunov dimension of 
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